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We analyze quantum algorithms for cloning of a quantum measurement. Our aim is to mimic two 
uses of a device performing an unknown von Neumann measurement with a single use of the device. 
When the unknown device has to be used before the bipartite state to be measured is available 
we talk about 1 — > 2 learning of the measurement, otherwise the task is called 1 — > 2 cloning of a 
measurement. We perform the optimization for both learning and cloning for arbitrary dimension 
of the Hilbert space. For 1 — > 2 cloning we also propose a simple quantum network that realizes the 
optimal strategy. 
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I. INTRODUCTION 



Arbitrary processing of a classical information can be 
described by strings of bits, and can be performed by 
a fixed device for example a processor of any personal 
computer. As a consequence, we do not need to build 
new devices for different computations, but we just need 
to copy bit strings carrying the appropriate program. 
Situation dramatically changes, when the systems carry- 
ing the information are governed by quantum mechanics. 
Unknown states of quantum systems can not be copied 
perfectly pQ and the no-programming theorem [2] pre- 
vents existence of universal quantum processors. This 
means that quantum programs can not be copied and 
that by using registers of qubits (two level quantum sys- 
tems) one can not deterministically realize all quantum 
information processing functions with a fixed processor. 
So in contrast to classical devices, quantum ones can- 
not be replicated by just copying the program for them. 
Copying of quantum states was extensively investigated 
[3HZ]. On the other hand copying of quantum devices 
did not receive so much attention even though it is a 
fundamental and equally important quantum informa- 
tion processing task. Similarly to states, quantum trans- 
formations are often used in quantum key distribution 
schemes [5HTT] to encode bits, so analysis of possible at- 
tacks by cloning them are needed. Cloning of transfor- 
mations was yet analyzed only for the case of unitary 
transformations [11] . In the present paper we investigate 
cloning of measurement devices, which can be seen as a 
cloning of certain measure-and-prepare transformations. 
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More precisely, when a measurement is an intermediate 
step of a quantum procedure its outcome can influence 
the following operations. This feed forward of the classi- 
cal outcome can be conveniently described using a quan- 
tum system into which the outcome is encoded into per- 
fectly distinguishable orthogonal states. In this sense a 
quantum measurement with only classical outcomes can 
be seen as a channel, which first measures the input sys- 
tem and based on the outcome prepares a state from a 
fixed orthogonal set. 

The term cloning of observables has been used in Ref. 
|12) referring to state cloning machines preserving the 
statistics of a class of observables. In the present paper 
the objective is to actually mimic two uses of an unknown 
measurement device, while using it only once. We would 
like to construct a replication strategy that would work 
for arbitrary von Neumann measurement E, even if it is 
provided as an unknown black box. The most natural 
operation of a replication strategy is based on modifying 
the bipartite state before E is actually used. In this case 
we talk about 1 — > 2 cloning of a measurement device. 
The most general representation of any cloning strategy 
is depicted below 




(the double wire carries the classical outcome of the mea- 
surement). 

On the other hand, one might ask how well the task can 
be accomplished when we use the measurement, before 
we have an access to the bipartite state of interest. We 
denote this scenario as 1 — >■ 2 learning, and any learning 
strategy can be depicted as follows 
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(2) 



In the present paper we will analyze only the above two 
scenarios, even though one can think of more general 
versions of the problem, where the M replicas have to 
be produced out of N uses of a measurement device. For 
example N — > 1 learning was analyzed in Ref . |13| . From 
comparison of Eqs. ([I]) and ^ one can see that learning 
is a particular instance of cloning in which the first step 
is restricted. That being so, it is clear that the perfor- 
mance of the optimal learning cannot be better than the 
performance of the optimal cloning. 

The paper is organized as follows. In Sec. [Il]we expose 
the formulation of the optimal learning and cloning in 
mathematical terms. In Sec lIIII we review the framework 
of quantum combs that is used as main tool throughout 



the paper. In Sec. IV the problem is simplified exploiting 
all the symmetries that can be useful. Sections [V] |VI| are 
devoted to derivation of optimal cloning and learning, 
respectively. The paper is closed by concluding remarks 
in Sec. ECD 



II. MATHEMATICAL FORMULATION OF THE 
PROBLEM 

Let us now formulate the problem mathematically. 
First of all, we should be able to evaluate the performance 
of the chosen replication strategy TZ. Hence, we need a 
quantity that expresses the closeness of a replicated mea- 
surement to a desired bipartite von Neumann measure- 
ment. In the following Lemma we introduce a function 
J^(P, Q) that quantifies the closeness of a POVM Q to a 
von Neumann POVM P. Throughout the paper we shall 
use the bold face notation for objects that are composed 
from several elements. For example P = {P}f =1 denotes 
the POVM with elements Pi and J = {/} is the single 
outcome POVM. 

Lemma 1 (Fidelity criterion for POVM) Let £ := 

{1, . . . , d} be a finite set of events and P C £(Jt?) and 
Q C £(Jt?) be two POVM's, such that one of them is a 
von Neumann measurement. Consider now the quantity 



d 

J?(P,Q) := - J] Tr [PiQi]. 



(3) 



i=l 



Then ^ = l^P i = Q i ^i and & < 1. 

Proof. Without loss of generality we can assume that 
P is a von Neumann measurement and that we have 



Pi — \i)(i\ where \i) is an orthonormal basis of Jff. Then 
for Qi = Pi = \i)(i\ we have 

d 1 d 

* = ^E^Q.] = rfETWH = !■ ( 4 ) 

i—1 2—1 



On the other hand if & = 


1 we have 
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(5) 



which implies Yli^tj (*l Qj K) = 0- Since Qj > 0, we 
must have (i \ Qj \i) = for all i ^ j, and consequently 
Qj = otj\j){j\ with ctj > 0. Finally the condition 

Yfj=i a j\j)U\ = Sy=i Qj = 1 implies a 3 = 1 and thus 
Qj = Pj . Proving that & < 1 is easy. Since Qi is an ele- 
ment of a POVM we have (i\ Qi \i) < 1 and consequently 

& = i Eti <t\ Qi N) < i- ■ 

Since we assume that the unknown measurement E to 
be replicated is a von Neumann POVM, we can write it 
in the following form 



Ei = \4n){4>i 



(6) 



where {\4>i)}f = i is an orthonormal basis of the Hilbert 
space W. All the POVM's of this kind can be generated 
by rotating a reference POVM {|i)(i|}f =1 by elements of 
the group of unitary transformations SU(rf) as follows 



Ef ] = U\i)(i\tf UeSU(d). 



(7) 



Let us denote the bipartite POVM replicated by the 
strategy U as G {u) = G{1Z,E {U) ). Our task is to find 
such replicating strategy "K. that the elements are 



as close as possible to E^ (g> E\ u> ■ Assuming that the 

unknown POVM E^ U ' is randomly drawn according to 
the Haar distribution, we choose the quantity: 

F[R]:= J dU ^(G (U) ,E^ U) (g> E {u) ) = (8) 



as a figure of merit for the replicating strategy. Hence, 
after choosing one of the two considered scenarios (1 — > 2 
cloning or learning) the goal is to find a strategy 1Z, that 
maximizes F\TV\- 
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III. PRELIMINARY CONCEPTS 

In this section we introduce the necessary notation 
and review the general theory of Quantum Networks, 
as developed in [T31 US]. Let us first recall the Choi- 
Jamiolkowsky isomorphism. It is an isomorphism con- 
necting any quantum operation (i.e. completely positive 
map) M. : B(J^f in ) H> B{,¥C out ) to a positive operator 
M e B{J^ ou t <S> fflin) defined as follows: 

M := M®l{\w){u\) (9) 

where T is the identical map on B(Jtfi n ), \lu) := 

S n l n ) \ n ) G ® ^™ ano - we fi xe d an orthonormal 
basis {\n)} on The action of M on a given input 

state p can be expressed in terms of M as: 

M{p) = Tr in [M(I®p T )\ (10) 

where Tr^„ denotes the partial trace over M[ n and the 
superscript T marks the transposition with respect to 
the basis {|n)}. 

Under the term Quantum Network we mean a network 
of quantum devices part of whose inputs and outputs are 
connected, while the remaining ones are forming open 
slots of the network into which sub-circuits can be later 
inserted. A network with (N — 1) open slots has N input 
and N output systems, that we label by even numbers 
from to 2N — 2 and by odd numbers from 1 to 2N — 1, 
respectively. Each network can be visualized as in Eq. 

1 2 3 2N- 2 2J V-1 
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where the wires represent the connections of output sys- 
tems to next inputs. This flow of quantum systems in- 
duces a causal order among the wires , according to which 
the input system m cannot influence the output system 
n if m > n. 

A Quantum Network TZ can be represented in terms 
of its Choi-Jamiolkowsky operator R, called quantum 
comb, which is a positive operator acting on the Hilbert 
space M'out ® Mi n where J^ ut := ®^Z<z ^2j+i, -ffin ■= 
®j=0 ^2j> an d •ffi'n being the Hilbert space of the n-th 
system. For a deterministic quantum network (i.e. a net- 
work of quantum channels) the causal structure implies 
the following normalization condition 

Ti 2k -i[R (k) } = h k -2®R {k - 1] k = l,...,N (12) 

where B.W = R, R^ = 1, i?« e -C(®^ Jtk), TWi 
denotes the partial trace on J^ik—i an d lik-i is an iden- 
tity operator on fflik-i- 

We can also consider probabilistic quantum networks 
(i.e. networks of quantum operations), whose Choi- 
Jamiolkowsky operators must satisfy 

< R < S (13) 



where S is the Choi-Jamiolkowsky operator of a deter- 
ministic network. 

We call generalized instrument a set of probabilistic 
quantum networks TZ := {TZi} such that the set R := 
{Ri} of the corresponding Choi operators satisfies 

^Ri = Rn (14) 

i 

where Rq is the Choi operator of a deterministic network. 

Two quantum networks TZi and 7?. 2 can be connected 
by linking some outputs of TZi (72-2) with inputs of TZ 2 
(TZi), thus forming a new network IZ3 := IZi * TZi. We 
adopt the convention that the wired to be connected are 
identified by the same label. The connection of the two 
quantum networks is mathematically represented by the 
link product of the corresponding Choi operators R± and 
i?2, which is defined as 

R 1 *R 2 =Tr J e[R e 1 *R 2 ], (15) 

9 jxr denoting partial transposition over the Hilbert space 
Jtf of the connected systems (recall that we identify the 
Hilbert spaces of connected systems with the same la- 
bels). 

As we pointed out in the introduction, the classical 
outcome of the inserted measurement can influence the 
next operation of the network. In order to take the feed 
forward of the classical outcome into account it is conve- 
nient to describe the measurement device to be replicated 
as a measure-and-prepare quantum channel 

d 

£( u \p) = Y / nEPp}\i)(i\, (16) 

which measures the POVM E^ u > on the input state and 
in the case of outcome i prepares the state \i) from a fixed 
orthonormal basis on the output of the channel. Within 
this framework the classical outcome is encoded into a 
quantum system by preparing it into a state from a set 
of orthogonal states. The Choi-Jamiolkowski representa- 
tion of the channel £^ u > is the following 

d t d 

e (u) = £ m\ ® Ef> = J2 m ® u*\i)(i\u T , 

i=l i=l 

(17) 

where X T denotes the transpose of X in the basis 

mytv 

Since we want the replicating network TZ to behave 
as two copies of the POVM E^ u ' upon insertion of a 
single use of £( u \ we have that TZ is actually a general- 
ized instrument R = {R%j}f J=1 where i,j is the couple 
of outcomes of the two replicated measurements. The 
normalization of the generalized instrument R = {Rij € 
C(JffA ® ® J^c ® J^h)} nas to obey the following 
equations: 
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1 — > 2 cloning 

^R ij = Rn = lD®S ABC Ty c [S] = Iab, (18) 

1 — > 2 learning 

Rij =Rn= Iabd <8 Pc Tr\p] = 1, (19) 



where the capital subscripts denote the Hilbert spaces 
on which the operators act and we use the labeling intro- 
duces in the Eqs. @, 

The replicated POVM is then equal to 



(20) 



£ (k\ c (k\ D (C/t ® J)^(C7 ® J) |fc) c 



IV. 



SYMMETRIES OF THE REPLICATING 
NETWORK 



which can be seen as FL 



Sij * M 



Sij * (see Eq. 



(17)) with the link performed on system D carrying the 
classical information. We can easily prove that 1Z is a 
generalized instrument. Indeed we have 



]T Ra = Y, S » * e{I) = Sn * , 



(24) 



where the link is performed only on the space J&h- The 



operator in Eq. ( 24 ) is the Choi- Jamiolkowski operator of 



a deterministic quantum network satisfying the same nor- 
malization conditions as Sq. Since we 
show that S and TZ produce the same replicated POVM 
G+P when linked with the single use of E^ u \ as follows 



13 * ^CD 



J CD 



P * E {U) 



((k\ D (k\ c U^)S ij (\k) D U\k) c ), (25) 



where the explicit form of the star product will be used 
later. The thesis then holds with R'^ k := (k\Sij\k)M 

The restriction to diagonal networks allows us to sim- 
plify the figure of merit (Eq. (fsl) ) as follows 



In this section we utilize the symmetries of the figure 
of merit ^ to simplify the optimization problem. These 
considerations apply both to cloning and learning of a 
measurement device. The first simplification relies on the 
fact that some wires of the network carry only classical in- 
formation, representing the outcome of the measurement. 
The classical information encoded in the choice of a state 
from basis {\i}} can be read without disturbance by the 
measure and prepare channel M with Choi- Jamiolkowski 
operator M = £W. Thus, inserting channel M. between 
the use of a measurement device E^ u ^ and the network 
TL will not change the operation of the scheme, i.e. 



-Te> m 



(21) 



As a consequence we have the following lemma. 



Lemma 2 (Restriction to diagonal network) The 

optimal generalized instrument 7Z, J2i j Rij = -^n 
maximizing Eq. fify can be chosen to satisfy: 



where < R' lj k € C{ J^ A ® J? B ® J^c)- 



(22) 



F[n] 



xU* m ®U\ijk) (ijk\ U T ® 2 <S> 



where \ijk) = \i) A \j) B \k) c and we applied Eq. (22). 

Since the performance of the scheme is evaluateaas an 
average over all possible "orientations" of the replicated 
measurement device, there exists a symmetrization pro- 
cedure that can make any strategy covariant (i.e. having 
property from Eq (27)), without affecting the figure of 
merit. 



U 



U 



u 



(27) 



This translates into mathematical terms as follows. 



Lemma 3 (Restriction to covariant networks) 

The operators R[j k that maximize Eq. (26) can be 
chosen to satisfy the commutation relation 



[R! ii<k ,U A ®U B ®Uc}=0. 



(28) 



Proof. Let Sij be the Choi representation of a general- 
ized instrument corresponding to a quantum network S. 
Let us define network 1Z as 



Rij :=J2(k\Sij\k)^\k){k\, 



(23) 



Proof. Suppose that the generalized instrument corre- 
sponding to k is optimal. Then one can easily check 
that also the instrument R'^ k defined as follows 



J& fc := / dU(U 



*®2 



WW^W) (29) 
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is suitably normalized and satisfies [R 1 ^ k ,U* ®U* ®U] — 
0. Generalized instrument TL corresponds to a strategy 
where random unitary U\ IP, U is applied before and 
after the original strategy S to systems A, B, C , respec- 
tively. From the integration in Eq. ( 26 1 it is obvious that 
the value of F for the above choice of R'^ k is the same 

as for S 'ij,k M 



The commutation relation ( 28 ) allows us to rewrite the 
figure of merit as 



F W = J» E W\ R kk \ijk) 



ABC ■ 



(30) 



i,j,k 



Another symmetry we can utilize is related to a simul- 
taneous relabeling of the outcomes of the inserted and 
produced measurements. We shall denote by a the ele- 
ment of §dt the group of permutations of d elements, and 
by T a the linear operator that permutes the elements 
of basis {\i}} according to this permutation, in formula 
T a \i) = \a(i)}. Let us note that the complex conjugation 
and transposition are defined with respect to the basis 
{!<)}, soT ff =T*. 

Lemma 4 (Relabeling symmetry) Without loss of 
generality we can assume that the operators R'^ k that 



maximize Eq. (26) satisfy the relation 



R', 



R' 



Hj,k L cr(ij,k) - 

where we shortened a (ij, k) := (o~(i) a(j),a(k)). 



(31) 



Proof. Suppose that network S characterized by op- 
erators Sij is optimal and satisfies both conditions (22) 
and ( 28 1 . Let us then define 



a£S d 



(32) 



a£S d 



where the last identity in (32) follows from the commu- 
tation relation (28) with U = T a , The operators 
correspond to a valid quantum network 7t, because H is 
a convex combination of networks Z a defined by Eq. ( 22 ) 
with Z'? k = S' a ^- k y Quantum network 1Z. operationally 
corresponds to relabeling of the outcomes of the inserted 
and replicated measurements by permutation a. The fig- 
ure of merit for 7?. is 



LEE <*(y*)i WW)) = m- 



d 2 d\ 



(33) 



cGSd i,j,k 



It is easy to prove that R'^ k satisfy Eq. (31 ) 



Remark 1 The properties (22), (28) and (31) induce the 



following structure of the replicated POVM's: 



G 



(t0 



(UT a 



(T^f 



(34) 



The advantage of using the relabeling symmetry is the 
reduction of the number of independent parameters of 
the quantum generalized instrument. Let us define the 
equivalence relation between strings ijk and i'j'k' as 



ijk r~j ijk ijk = <x(i j k ), 



(35) 



for some permutation a. Thanks to Eq. (31) there 
are only as many independent R\ - k as there are equiv- 
alence classes among sequences (ij,k). There are four 
or five equivalence classes depending on the dimension 
d being two or greater than two, respectively. We 
denote the set of these equivalence classes by L := 
{xxx, xxy, xyx, xyy, xyz}. 

Based on lemma|4]we can write the optimal generalized 
instrument as follows 



Rab 



R' 



ij,k 



R' 



a{ij,k) 



(36) 



where (ab, c) is a string of indices that represents one 
equivalence class from L. 

The figure of merit can finally be written as follows 



F[Tt\ = ^ n(ab,c)(R ab , c 



d 2 



(37) 



(ab.c)SL 



where n(ab, c) is the cardinality of the equivalence class 
denoted by (ab, c), and (R a b,c) = (ijk\ R'^ k \ijk) for any 
string ijk in the equivalence class denoted by (ab, c). As 
a consequence of Schur's lemmas, Eq. (28) implies the 
following structure for the operators R a b,c 



Rab,c = ®P V 



r ab,ci 



(38) 



where v labels the irreducible representations in the 
Clebsch-Gordan series of U* A ® Ug ® Uc, and P u acts 
as the identity on the invariant subspaces of the repre- 
sentations v, while r v ah acts on the multiplicity space of 
the same representation. 

Depending on the dimension d — 2 or d > 2 we have 
two different decompositions. In the former case, we have 



R a b,c — P a 



r ab,c + ^ r ffc,c' 



(39) 



where r" h c is a positive 2x2 matrix, while r ah c is a non- 
negative real number. The projections P^ on the invari- 
ant spaces of the representation U* ® U* <g) U are the 
following 



m—l 

pp = i®p+ - p a ®\+){+\, 



(40) 
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where W± = (\uj) \m)±\m) |u;))/[2(d±l)]5, and P+, P~, 
are the projections onto the symmetric and antisymmet- 
ric subspace, respectively. When d > 2, on the other 
hand, we have 



R ab ,c = P a ® r%_ c + P^ h r + prrZ, „, (41) 



ab,c ab,c^ 

Y*< rind nrSY 



where r° b is a positive 2x2 matrix, while r^ b and r , 

are non-negative real numbers. The projections P^ on 
the invariant subspaces are the following 

d 

i~ ® |a)<6| = £ |*° ))<<*U a, &€{+,-} 



P? = I®P+ -P a ®\+){+\, 
p 7 = i®p- -P a <g>|-)(-|. 



(42) 



The last symmetry we are going to introduce relies 
on the possibility to exchange the inputs (Hilbert spaces 
J^A and M'b) of the two replicated measurements with 
simultaneously exchanging their measurement outcomes, 
while the figure of merit is left unchanged. 



Lemma 5 The operators R ab ,c in Eq. (38) can be chosen 
to satisfy 

R abtC = SR ba , c S V(a6, c) G L (43) 
where S is the swap operator S \k) A \j) B = \j) A W) b- 

Proof. The proof can be done by the following averag- 
ing argument. Let us define Rij t k '■= \{R'ij k + *>i2^ fe S). 
It is easy to prove that {Rjj ,k} satisfies the correspond- 
ing normalization (Eq. ( 18 ) for cloning or Eq. ( 19 ) for 



learning) and that it gives the same value of F as R[ 



ij,k- 



Eq. (W3| together with the decomposition ( 41 ) gives 
for y(ab~c) € L 



o,r nh „cr 2 = r h 



P 7 7 



(44) 



where cr 2 



1 

-1 



As a consequence of Eq. ( 38 1 , the figure of merit in 
Eq. ( 37 ) can be written as 



F =4 E n(ab,c)Tr[\ijk){ijk\J2P V ®<b,c] 



where 



d 2 

(ah,c)GL 

E \ E ^\K b Ac] 

v (ab,c)el 

--F a + F +F 1 



n{ab,c) v 



(45) 



(46) 
(47) 



and ij, k is any triple of indices in the class denoted by 
ab,c. Notice that n(xx,x) — d, n(xx,y) — n(xy,x) = 
n(xy, y) = d(d — 1), n(xy, z) = d(d — l)(d — 2) and in the 
case d — 2 F 1 — (i.e. does not appear). 
In particular, by direct calculation we have 



A Q = [ d +! ° ) A 

— Vr.r..r Q 



xx,y 



d+l 








a 

xx,y^z 




d 




2(d + 


1) 



A' 3 =— , A^ =A' 3 

xx, x ^ _|_ 1 xx,y xy,x 

A' 3 =1 A' 3 =- 

xy,y i xy,z ^ 5 

A 7 = A 7 =0 A 7 = A 7 = d ~ 2 

xx,x xy,y > xx,y xx,y 2(ti 1) 

A 7 = l - 

xy,z 2 



OPTIMAL CLONING 



(48) 



In this section we turn our attention to the cloning sce- 
nario. Cloning is less restrictive than learning, since we 
allow the two states to be measured to be available at the 
same time as the single use of the measurement device. 
The normalization condition for the 1—^2 cloning reads 



ij,k 



\k}(k\ D R' i:j k = I D ® Sabc Tr c [S] = I AB 



(49) 



which implies the following 

I AB =Tr c [R xx , x ] + (d-l){d-2)Tr c [R xy , z ] 

+ (d- l)Tlc[Rxx,y + Rxy,x ~t~ Rxy,y\- (^0) 

From the commutation [R ab , c , U A ®U B ®Uc] it follows 
that [Trc[Rab,c], U A ® U* B ] = and taking the decompo- 
sition R ab , c = Y. v pv 



the normalization constraint ( 50 1 becomes 



P 



Trn\P" 



=^E z. »<■ 

v (ab,c)eL 



ab c along with definition (47 1, 
50 1 becomes 

]P ± . 



^ab,c\ 



(51) 



We take the trace of the previous equation to obtain 
the following equivalent formulation of the normalization 
constraints 



d+ = d a ^ s ab ,c 



b ab,c> 



(ab,c)£L 



da S ~ 
(ab,c)eL 

d 1 S 



7 



;6,c> 



(52) 
(53) 



d -=d a }^ s o 

(a6,c)GL (ab.c)eL 

where d± = Tr[P±], d v = Tr[P u ]. If we introduce the 
notation 

J 

'a.l 

0, 



J .- I S 'a,bc 



/0 

si i. := I n ,7 

a.bc 



s 



n + = 







n - 



o o 
o 1 



(54) 



7 



the normalization constraints ( 52 1 and ( 53 ) can be rewrit- 
ten as 



i i/,(a,bc)£l_ 



IT 



d+ 




n 



E 



'(a, he) 



IT 



v v,{a,bc)e\- 





d_ 



(55) 



In order to solve the optimization problem we have to find 
the set s := {s%,£ £ L,i/£ s» E C(C 2 ),s^ > 

subjected to the constraint ( 55 1 that maximizes the figure 



of merit (45); we will denote as M the set of all the s 



satisfying Eq. (55 1. Since the figure of merit (45) is 



linear and the set M is convex, a trivial result of convex 
analysis states that the maximum of a convex function 
over a convex set is achieved at an extremal point of the 
convex set. We now give two necessary conditions for a 
given s to be an extremal point of M. Let us start with 
the following 

Definition 1 (Perturbation) Let s be an element of 
M. A set of hermitian operators z := {z^} is a per- 
turbation of s if there exists e > such that 



hz e 



V/iG 



(56) 



where we defined s + hz := {s^ + hz%\h £ [— e, e]}. 



By the definition of perturbation it is easy to prove that 
an element s of M is extremal if and only if it admits 
only the trivial perturbation Zg — W, v. We now ex- 
ploit this definition to prove two necessary conditions for 
cxtrcmality. 

Lemma 6 Let s be an extremal element of M . Then Sg 
has to be rank one for all £, v . 



Proof. Suppose that there is a s v v 



{0,...,0,4/,0,...,0}, zi = 



€ S 



is an admissible 



a b 

"" ~ I c d 

which is not rank one; then there exist e such that z := 

1 
x 1 
perturbation. ■ 

The above lemma tells us that without lost of general- 
ity we can assume the optimal s to be a set of rank one 
matrices. Let us now consider a set s such that s\ is rank 
one for all t, v\ any admissible perturbation z of s must 
satisfy 



(57) 



(58) 



where the constraint (57) is required in order to have 



the normalization ( 55 ) . Let us now consider the map 



f-C{€ 2 ) 



f(A) - ( Tr[Il+A] 
JKA) — I Tr[n~^4] 



exploiting this definition Eq. (|58|) becomes 









(59) 



Suppose now that the set s has N > 3 non-zero elements; 
then {/(s£)} is a set of N > 3 vectors of C 2 that cannot 
be linearly independent. That being so, there exists a set 
of coefficients {c^} such that X^£ c f f( s e ) = and then 
Zp = c v l s v t is a perturbation of s. We have then proved 
the following lemma. 

Lemma 7 Let s be an extremal element of M. Then s 
cannot have more than 2 non-zero elements. 

Lemma [6] and Lemma [7] provide two necessary conditions 
for extremality that allow us to restrict the search of the 
optimal s among the ones that satisfy 

s = { si , si!, } Rnk(4' ) = Rnk(4" ) = 1 



IT 



IT 



(60) 



The set of the above s is small enough to allow us to 
compute the value of F for all the possible cases. It turns 
out that there are two choices achieving the highest value 
of fidelity 



4 

3d' 



(61) 



They are defined by s = {s% x x , s" y x } and s 
s £»,i/}' where 



- , ~~ 93~~ ] = A 

\ Q 



B = 




s a = B 



s xy,y = 0- z Ba z . 



(62) 



From the linearity of the link product and our figure of 
merit it follows that also any convex combination of the 
above two strategies will give the optimal performance. 
In the rest of the paper we consider the equal convex 
combination of the above two strategies 



,a -A 

°xx,x ' 



1 



:B, 



1 



-a z Ba z 



(63) 



xy.x 2 ' "ViV 2 

because it treats the two clones in the same way. Us- 



s v e + hz v t > 0, while Eq. ( 58 ) tells us that s + hz satisfies ing Eq. ( 25 1 one can derive the form of the replicated 



POVM corresponding to the above choice of the optimal 
generalized instrument. 



G 



1 



P+(Ef ) ®I B )P+ 



1 



9d(d+ 1) 

6', , = -^-[Q+(Pf } ® 7 S )Q+ + Q~{E) U > ® J B )Q-], 



where Q ± = 1/V9d(d+ 1) P+ ± l/\/2 P" 



A. Realization scheme for the optimal cloning 
Network 



In this section we describe the inner structure of the 
optimal cloning network. First we notice that the choice 
from Eq. ( 63 1 corresponds to the generalized instrument 



Pj 
Pj 



l*> <*l ® 



9cL 



2(d- 1) 



2(d-l) 

1 i _l\ /d_ i \ 



(64) 




9d 



The generalized instrument P can be realized by the 
following network 




(65) 



The first step consists of a control SWAP gate, which is 
described by the unitary 



Ucs — Ta^c 

+ / | . K 



1 T B ->k |u;\u|i- 



with the control qubit prepared in the state |+) = 
i(|0) + |1)). We defined T X ^ Y = Eil*)r<*lx and 
named Ml the 2-dimensonal Hilbert space of the con- 
trol qubit with {|0) , |1)} being an orthonormal basis on 

In the second step we have three commuting actions: 

• the single use of the measurement device E^ u ' is 
applied on system C and its outcome is recorded 
on a classical memory D 

• system K is discarded 



system L undergoes a 3-outcome measurement de- 
scribed by the POVM P defined as follows: 



Pi = 
P 2 = 



(9d(d + l)-2) 



-x- 



9d(d + l) 

XVI ft = owhMVK 



i 



9d(d+l) 



(66) 



The last step is just a classical processing / of the out- 
come k of the measurement E^ u ' and of the outcome n 
of POVM P. The function / that produces the outcome 
— f(k,n) of the network is defined as follows: 



f(k,n) 




(67) 



where the outcome j in the second and third case is ran- 
domly generated with flat distribution. 

In order to prove that this network is described by the 
generalized instrument in Eq ( 64 ) we first realize that the 
action of the POVM P and of the processing / can be 
represented by the bipartite POVM Q on systems D and 
L defined as 



(9d(d+l)-2) I 
9d(d+l) I 



if z 



J 



er z \ifi)(i/)\cr z 



d-1 



if i ^ j 



C , s D 



Q 



c 



L 



D 



Finally, one can check the identity 



Pj 



\+)(+\*\Ucs))((Ucs\*{Q 



I 



K) 



(68) 



(69) 



It is worth to notice that the optimal cloning of mea- 
surement device has some features in common with the 
optimal cloning of unitaries. Both in the cloning of uni- 
taries and in the cloning of von Neumann measurements 
the first step is to perform a control-SWAP of the two 
input systems with the control qubit prepared in the su- 
perposition ^/j(|0) + |1)). We could give an intuitive 
explanation of this feature in terms of quantum paral- 
lelism: for a bipartite input \x) the unknown mea- 
surement acts on both input states via a superposition 

%(\x) 2 \0 A + \0 2 \x) A )- 



VI. OPTIMAL LEARNING 

Our goal in this scenario is to create two replicas of the 
measurement after it was used once. Let us consider the 
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normalization constraint for the generalized instrument 
Rij. Since ^ 4 ■ Rij has to be a deterministic network, 
we have 

]T \k) {k\ D ® k = Iabd ® PC, Tr[p] = 1 (70) 

where p has to be positive operator. The commutation 
relation (28 1 implies [p, U] = and so we have p = \lc- 
Writing J ASCD as £ fc (A;| D _ ® (I mot ® P a + pP + P~<) 
we can rewrite the normalization conditions as follows 



(ah,c)£L 
(ah,c)£L 



1 . 



f = 0,7 



(71) 



Let us now maximize the figure of merit under these con- 
straints. The maximization of Fp and F 7 is simple and 
yields 



(72) 







d? 


P — 

7 2d 2 








X 


^xy : x 


- - 

J xy,y 


s' 3 


= 


s 7 

' J xx 


X 


- S 7 
xx,y 


_ „7 — 
u xy,x 


s 7 

a xy,y 


= 


' J xx 


y 


- s 7 

X V: Z 


1 

~~ d' 







(73) 



Let us now consider the maximization of F a ; the normal- 
ization constraint for the a subspace gives 



/ j *ab,c ~~ rfi a ab,c ~ u ' 

(ab.c)eL 

E c: c ' + = o. (74) 

(ab,c)£L 



(a&,c)eL 

Ea, 
S ab 

(afc,c)GL 



d' 
1 



Inserting the explicit expression of the A" b into Eq. ( 45 ) 



and taking into account Eq.(44) we have 



dF„ =Tr 



Tr 



■ s a,+,+ s a, + ,— 



■ „a, + ,+ „a, + ,— 

°xy,x °xy,x 

„a,— ,+ „«,—,— 

• °xy,x °xy,x 




2d(d 2 - 1) d + 1 



where in the derivation of the bound (75) we used the 



positivity of s" x and the constraints (74 1. The upper 



bound ( 75 1 can be achieved by taking 



J xy,z 



xx,y 



-2a 






xy.x 



2d" 2d 



(76) 



where wc defined a := s"^ ,+ . Eq. (75) gives the value 
of F a as a function of a; the maximization of F a (a) with 
the constraint < a < ^ is easy and gives 



F a = 



4(2d- 1) 



for a = 



1 



3d 2 (d 2 -l) ~ 18d(d-l)' 

and then for d > 3 we have 

9d 2 + 16d - 17 



6d 2 (d 2 - 1) 



3 

2d 2 ' 



(77) 



(78) 



For d — 2 the invariant subspace ML does not appear 
and the fidelity becomes F = F a + Fp = X. 

Using Eq. (25) it is possible to derive the form of the 



replicated POVM corresponding to the optimal general- 
ized instrument. 



16d-2 



9d(d 2 
1 

d- 1 



1 



P+(K 



(u) 



q'+(eP 



Ih)P^ 



d(d 2 - 1) 



P' 



Ib)Q' + + Q'-{Ef ] ®I B )Q'- 



d(d-l) 2 (d-2) 
d-3 



E 



(U) 



Ib)P~ 



(d-l) 2 (d-2) 



where Q ,zk = l/^/9d(d- 1) (P+ ± 3 P~). 

One can now compare the performance of the optimal 
1 — s- 2 cloning and learning. The optimal values of F 
depending on the dimension d are plotted on Figure [T] 
As expected the optimal cloning strategy largely outper- 
forms the optimal learning strategy with a fidelity, which 
is a factor d larger, as one can see from Eqs. (61) and 
(78). Similar distinction arises also for comparison of 



cloning and learning of unitary channels (for details see 
refs [H]). It is also worth noting that the optimal learn- 
ing strategy achieves a greater fidelity than the incoher- 
ent strategy in which one first make the optimal estima- 
tion of the measurement and then conditionally prepares 
two copy of the estimated measurement (one can prove 



that for this last strategy one has F„ 



— ( d+2 
\d(d+l) 



VII. CONCLUSIONS 



In the present paper we focused on 1 — > 2 cloning and 
1 — > 2 learning of von Neumann measurements. Even 
though both problems can be easily formulated in the 
usual language of quantum mechanics, the necessity to 
handle the measurement outcome in the remaining part 
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FIG. 1: Optimal 1 — > 2 cloning and learning of a measure- 
ment device: we present the values of F for different values 
of the dimension d. The squared dots represent the optimal 
1 — > 2 cloning, the round dots represent the optimal 1 — > 2 
learning and lowest curve corresponds to a strategy in which 
one performs the optimal estimation followed by the prepara- 
tion of the estimated measurement. 



sure&prepare channel and we employed framework of 
quantum combs to perform the network optimization. 
Thanks to symmetries of the figure of merit the prob- 
lem was simplified and solved for arbitrary dimension 
of the measurement's Hilbert space d. In section (VA) 
we proposed a realization of optimal 1 — > 2 cloning of 
measurements. The proposed scheme has some similari- 
ties to optimal cloning of unitary transformations, since 
they both begin by the control-swap operation, which re- 
flects the presence of quantum parallelism. In this paper 
we expolit the measure&prepare representation of von 
Neumann measurement that allowed us to deal with feed 
forward of classical information in quantum networks. 
These tool could be in principle used to tackle other 
quantum information processing tasks in which classi- 
cal information is involved e.g. estimation and cloning of 
quantum instruments. 
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